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Abstract 


Two searches for physics beyond the standard model in events containing photons are 
presented. The data sample used corresponds to an integrated luminosity of 19.7 fb 1 
of proton-proton collisions at y/s = 8 TeV, collected with the CMS experiment at the 
CERN LHC. The analyses pursue different inclusive search strategies. One analysis 
requires at least one photon, at least two jets, and a large amount of transverse mo¬ 
mentum imbalance, while the other selects events with at least two photons and at 
least one jet, and uses the razor variables to search for signal events. The background 
expected from standard model processes is evaluated mainly from data. The results 
are interpreted in the context of general gauge-mediated supersymmetry, with the 
next-to-lightest supersymmetric particle either a bino- or wino-like neutralino, and 
within simplified model scenarios. Upper limits at the 95% confidence level are ob¬ 
tained for cross sections as functions of the masses of the intermediate supersymmet¬ 
ric particles. 
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1 Introduction 

Supersymmetry (SUSY) JTHZl is a popular extension of the standard model, which offers a so¬ 
lution to the hierarchy problem (SJ by introducing a supersymmetric partner for each standard 
model particle. In models with conserved R-parity |9j[T0|, as are considered here, SUSY parti¬ 
cles are produced in pairs and the lightest supersymmetric particle (LSP) is stable. If the LSP 
is weakly interacting, it escapes without detection, resulting in events with an imbalance p™ ss 
in transverse momentum. Models of SUSY with gauge-mediated symmetry breaking lHTI - fl7l 
predict that the gravitino (G) is the LSP. If the next-to-lightest SUSY particle is a neutralino 
(Yi) with a bino or wino component, photons with large transverse momenta (fj) may be pro¬ 
duced in —>■ yG decays. The event contains jets if the y ( | originates from the cascade decay 
of a strongly coupled SUSY particle (a squark or a gluino). 

In this paper, we present two searches for gauge-mediated SUSY particles in proton-proton 
(pp) collisions: a search for events with at least one isolated high-px photon and at least two 
jets, and a search for events with at least two isolated high-px photons and at least one jet. The 
discriminating variables are E™ 1SS for the single-photon analysis, and the razor variables Mr 
and R 2 [18, 19j| for the double-photon analysis, where Ey ,ss is the magnitude of p| 11,ss . These 
studies are based on a sample of pp collision events collected with the CMS experiment at the 
CERN LHC at a center-of-mass energy of 8 TeV. The integrated luminosity of the data sample 
is 19.7 fb _1 . 

Searches for new physics with similar signatures were previously reported by the ATLAS and 
CMS collaborations at y/s = 7TeV, using samples of data no larger than around 5 fb 1 l2Q] - j23| . 
No evidence for a signal was found, and models with production cross sections larger than 
~ 10 ftU 1 were excluded in the context of general gauge-mediation (GGM) SUSY scenarios [24]- 

m. 

This paper is organized as follows. In Section[2]we describe the CMS detector, in Section[3]the 
benchmark signal models, and in Section [4] the part of the event reconstruction strategy that is 
common to the two analyses. The specific aspects of the single- and double-photon searches 
are discussed in Sections [5] and [6j respectively. The results of the analyses are presented in 
Section [8j A summary is given in Section [9] 

2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam¬ 
eter, providing a magnetic field of 3.8 T. Within the superconducting solenoid volume are a 
silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and 
a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap 
sections. Muons are measured in gas-ionization detectors embedded in the steel flux-return 
yoke outside the solenoid. Extensive forward calorimetry complements the coverage provided 
by the barrel and endcap detectors. 

Events are recorded using a trigger that requires the presence of at least one high-energy pho¬ 
ton. This trigger is utilized both for the selection of signal events, and for the selection of 
control samples used for the background determination. The specific trigger requirements for 
the two analyses are described below. Corrections are applied to account for trigger inefficien¬ 
cies, which are evaluated using samples of data collected with orthogonal trigger conditions. 
A more detailed description of the CMS detector, together with a definition of the coordinate 
system used and the relevant kinematic variables, can be found in Ref. Il30l . 



2 


3 SUSY benchmark models 


3 SUSY benchmark models 

The two searches are interpreted in the context of GGM SUSY scenarios f24T(29fl . and in terms 
of simplified model spectra (SMS) scenarios !(3lTl34l inspired by GGM models. In these sce¬ 
narios, R-parity is conserved and the LSP is a gravitino with negligible mass. Four models are 
considered: 

GGMbino model In this model, squarks (q) and gluinos (g) are produced and decay to a final 
state with jets and a bino-like y\’. This production process dominates over electroweak 
production in the squark- and gluino-mass region accessible to the analyses. The mass 
is set to 375 GeV, leading to a Xi ~> Gy branching fraction of about 80% l26l . The events 
are examined as a function of the squark and gluino masses. All other SUSY particles 
have masses set to 5 TeV, which renders them too heavy to participate in the interactions. 
In most cases, the final state contains two photons, jets, and E™ 1SS . 

GGMwino model This model is similar to the GGMbino model, except that it contains mass- 
degenerate wino-like yi and yy particles instead of a bino-like x^. The common mass of 
the x i an d X \ is se t to 375 GeV. The final state contains a 77 , 7 V, or VV combination in 
addition to jets and E™ ISS , where V is a Z or W boson. With a 7 ) —> G 7 branching fraction 
of about 28%, approximately 48% of all events contain at least one photon. 

T5gg model This SMS model is based on gluino pair production, with the gluinos undergoing 
a three-body decay to qqyj, followed by y^ —> G 7 . All decays occur with a branching 
fraction of 100%. The final state contains at least two photons, jets, and Ej ISS . 

T5wg model This SMS model is also based on gluino pair production, with one gluino under¬ 
going a three-body decay to qqylj 1 , followed byy'j 1 -+ G 7 , and the other gluino undergoing 
a three-body decay to qqy^, followed by yj —» GW ± . All decays occur with a branching 
fraction of 100%. The final state contains at least one photon, jets, and E™ ISS . 

Typical Feynman diagrams corresponding to these processes are shown in Fig. |l] Note that 
for the two GGM models, the events can proceed through the production of gluino-gluino, 
gluino-squark, or squark-squark pairs. 

Signal events for the GGM models are simulated using the PYTHIA 6 Il35l event generator. 
The squark and gluino masses are varied between 400 and 2000 GeV. Eight mass-degenerate 
squarks of different flavor (u, d, s, and c) and chirality (left and right) are considered. The 
production cross sections are normalized to next-to-leading order (NLO) in quantum chromo¬ 
dynamics, determined using the PROSPINO Il36l program, and is dominated by gluino-gluino, 
gluino-squark, and squark-squark production. 

The SMS signal events are simulated with the MadGraph 5 If37j Monte Carlo (MC) event 
generator in association with up to two additional partons. The decays of SUSY particles, 
the parton showers, and the hadronization of partons, are described using the PYTHIA 6 pro¬ 
gram. Matching of the parton shower with the MADGRAPH 5 matrix element calculation is per¬ 
formed using the MLM 1381 procedure. The gluino pair-production cross section is described 
to NLO+NLL accuracy It36l(39i - (42'l . where NLL refers to next-to-leading-logarithm calculations. 
All SUSY particles except the gluino, squark, LSP, and Xi states are assumed to be too heavy 
to participate in the interactions. The NLO+NLL cross section and the associated theoretical 
uncertainty Ii43l are taken as a reference to derive exclusion limits on SUSY particle masses. 
Gluino masses of 400 (800) to 1600 GeV, and Xi masses up to 1575 GeV, are probed in the T5wg 
(T5gg) model. 
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Figure 1: Typical Feynman diagrams for the general gauge-mediation model with bino- (top 
left) and wino-like (top right) neutralino mixing scenarios. Here, the can decay to G 7 or GZ, 
with the branching fraction dependent on the x ® mass. The diagrams for the T5gg (bottom left) 
and T5wg (bottom right) simplified model spectra are also shown. 

For all the signal models, detector effects are simulated through a fast simulation of the CMS 
experiment f44l . 


4 Event reconstruction 

The events selected in this study are required to have at least one high quality reconstructed 
interaction vertex. The primary vertex is defined as the one with the highest sum of the pj 
values of the associated tracks. A set of detector- and beam-related noise cleaning algorithms 
is applied to remove events with spurious signals, which can mimic signal events with high 
energetic particles or large E™ 1SS 1451 1461. 

Events are reconstructed using the particle-flow algorithm B7ll48ll . which combines informa¬ 
tion from various detector components to identify all particles in the event. Individual particles 
are reconstructed and classified in five categories: muons, electrons, photons, charged hadrons, 
and neutral hadrons. All neutral particles, and charged particles with a track pointing to the 
primary vertex, are clustered into jets using the anti-kj clustering algorithm [|49l , as imple¬ 
mented in the FASTjET package B50l . with distance parameter of 0.5. The momenta of the jets 
are corrected for the response of the detector and for the effects of multiple interactions in the 
same bunch crossing (pileup) ll5ll . Jets are required to satisfy loose quality criteria that remove 
candidates caused by detector noise. 

Photons are reconstructed from clusters of energy in the ECAL H52j . The lateral distribution of 
the cluster energy is required to be consistent with that expected from a photon, and the energy 
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detected in the HCAL behind the photon shower cannot exceed 5% of the ECAL cluster energy. 
A veto is applied to photon candidates that match hit patterns consistent with a track in the 
pixel detector (pixel seeds), to reduce spurious photon candidates originating from electrons. 
Spurious photon candidates originating from quark and gluon jets are suppressed by requiring 
each photon candidate to be isolated from other reconstructed particles. In a cone of radius 
A R = V (A ;/) 2 + (At/j ) 2 = 0.3 around the candidate's direction, the scalar pj sums of charged 
hadrons (I n ), neutral hadrons (/ n ), and other electromagnetic objects (J 7 ) are separately formed, 
excluding the contribution from the candidate itself. Each momentum sum is corrected for the 
pileup contribution, computed for each event from the estimated energy density in the (t],(p) 
plane. Selected photons are required to be isolated according to criteria imposed on I n , I n , and 
ly as defined in Ref. H52l . 

5 Single-photon search 

The single-photon analysis is based on a trigger requiring the presence of at least one pho¬ 
ton candidate with pj > 70 GeV. The trigger also requires Hj > 400 GeV, where Hj is the 
scalar sum of jet pj values for jets with pj > 40GeV and \rj\ < 3, including photons that are 
misreconstructed as jets. 

In the subsequent analysis, we make use of the variable pj, which is defined by considering 
the photon candidate and nearby reconstructed particles, clustered as a jet as described in Sec¬ 
tion |4j If a jet (possibly including the photon) is reconstructed within A R < 0.2 of the photon 
candidate and the pj value of the jet is less than three times of that of the photon candidate 
itself, it is referred to as the "photon jet". If such a jet is found, p 7 is defined as the pj value 
of the photon jet. Otherwise, pj is the pj value of the photon candidate. We require photon 
candidates to satisfy pj > 110GeV and \p\ < 1.44. Also, in the subsequent analysis, we make 
use of the variable Hj, defined as for Hj in the previous paragraph but including the Pj val¬ 
ues of all selected photon candidates. The variables pj and Hj reduce differences between 
photon candidates selected with different isolation requirements compared to the unmodified 
variables pj and Hj. We require events to satisfy Hj > 500 GeV. The sample of events with 
isolated photons so selected is referred to as the 7ti g ht sample. The trigger efficiency for the 
selected events to enter the sample is determined to be 97%, independent of p 7 and Hj. 

We require at least two jets with pj > 30GeV and \rj\ < 2.5. The jets must be separated by 
A R > 0.3 from all photon candidates, to prevent double counting. In addition, the requirement 
Ej' ss > 100GeV is imposed and events with isolated electrons or isolated muons are vetoed. 
The selection is summarized in Table [T] Note that 0.16% of the selected events contain more 
than one photon candidate. 

The relevant sources of background to the single-photon search are: 

• multijet events with large E™ ISS , originating from the mismeasured momenta of some 
of the reconstructed jets. This class of events contains both genuine photons and 
spurious photon candidates from jets. This is by far the largest contribution to the 
background. 

• events with genuine Ej ,ss originating from the leptonic decay of W bosons, both 
directly produced and originating from top quark decays, which we refer to as elec- 
troweak (EW) background. 

• rare processes with initial- or final-state photon radiation (ISR/FSR), such as 7 W, 

7 Z (especially 7 Z —> yvv), and 7 ft production. 
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Table 1: Summary of the single-photon analysis selection criteria. 


Selection criteria 

T tight 

signal reg. 

Tloose 

control reg. 

Tpixel 

control reg. 

Isolation requirement 
Pixel seed 

Tight 

Vetoed 

Loose 

Vetoed 

Tight 

Required 

Trigger 

y-Hj trigger with 
p 7 > 70 GeV, H T > 400 GeV 
(using > 40 GeV, \ t]\ < 3.0) 

Photon(s) 

>1, p * 7 > 110GeV, \rj\ < 1.44 

Jet(s) 

>2, p J P ts1,2 > 30 GeV, |j/| < 2.5 

TT* 

il'p 

> 500GeV 

(using p ^ ts ' 7 > 40 GeV, 7/ < 3.0) 

Isolated e,// 

veto, pj > 15GeV, |// e 00 

< 2.5 (2.4) 

£miss 

E™ 1SS > 100 GeV (six ranges in E™ 1SS ) 


The kinematic properties of the multijet background are estimated from a control sample of 
photon candidates with isolation-variable values (l n , I n , h f ) too large to satisfy the signal pho¬ 
ton selection. We refer to these events as the 7i 00 se sample. Photon candidates of this kind 
typically originate from jets with anomalous fractions of energy deposited in the ECAL. Other 
than the orthogonal requirement of a 7 i 00 se rather than a 7 ti g ht candidate, events in this con¬ 
trol sample are selected with the same requirements as the 7 ti g ht sample, as summarized in 
Table [T[ Despite the different isolation requirement, this sample has properties similar to those 
of the 7 ti g ht sample, due to the use of p'j rather than photon pj in the definition of the event 
kinematic variables. Moreover, events in the 7i 00 se control sample are corrected for a residual 
difference with respect to the 7 t i g ht sample in the distributions of pj and hadronic recoil pj, 
estimated from events with E™ ISS < 100 GeV. The corrected distribution of a given kinematic 
property (e.g., E™ 1SS ) for 7 i oose events provides an estimate of the corresponding distribution 
for 7 t i g ht events. The uncertainty in the correction factors, propagated to the prediction, is fully 
correlated among bins in the signal region and is treated as a systematic uncertainty in the 
background yield. The limited statistical precision of the control sample dominates the total 
systematic uncertainty. Figure 2] (left) shows the E™ 1SS distribution from the 7a g ht sample and 
the corresponding prediction from the 7 i 00 se sample, for simulated multijet and 7 + jet events. 
No discrepancy is observed within the quoted uncertainties. 

The EW background is characterized by the presence of an electron misidentified as a pho¬ 
ton. The kinematic properties of this background are evaluated from a second control sample, 
denoted the 7 p i xe i sample, defined by requiring at least one pixel seed matching the photon 
candidate but otherwise using the 7ti g ht selection criteria, as summarized in Table [lj Events in 
the 7 P ixel sample are weighted by the probability / e -> 7 for an electron to be misidentified as a 
photon, which is measured as a function of the 7 candidate pj, the number of tracks associated 
with the primary vertex, and the number of reconstructed vertices in an event by determining 
the rate of events with reconstructed e 7 p i xe i and e 7 t j g ht combinations in a sample of Z —> e 1 e 
events. The event-by-event misidentification rate is about 1.5%, with a weak dependence on 
the number of vertices. A systematic uncertainty of 11% is assigned to / c _ >7 to account for 
the uncertainty in the shape of the function and for differences between the control sample 
in which the misidentification rate is calculated and the control sample to which it is applied. 
The predicted £ p ISS distribution for the EW background, obtained from a simulated sample of 
W boson and tt events, is shown in Fig. [ 2 ] (right) in comparison with the results from the direct 
simulation of events with 7 t i g ht originating from electrons. The distributions agree within the 
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19.7 fb 1 (8 TeV) 19.7 fb' 1 (8 TeV) 



Figure 2: Tests of the background estimation method for the single-photon analysis using simu¬ 
lated events in bins of E™ ISS . The direct simulation of 7 t i g ht events is compared to the prediction 
of the multijet background from simulated 7i 00 se events (left). Simulated events with 7ti g ht 
originating from generated electrons are compared to the simulated prediction using the EW 
background method (right). The blue hatched area represents the total uncertainty and is the 
quadratic sum of the statistical (red vertical bars) and systematic (red hatched area) uncertain¬ 
ties. In the bottom panels, the ratio of the direct simulation to the prediction is shown. 

quoted uncertainties. 

The contribution of ISR/FSR background events is estimated from simulation using leading- 
order results from the MadGraph 5 MC event generator with up to two additional partons, 
scaled by a factor of 1.50 ± 0.75 including NLO corrections determined with the MCFM H53l l54j| 
program. 

The measured E™ ISS spectrum in the 7ti g ht sample is shown in Fig. |3j in comparison with the 
predicted standard model background. A SUSY signal would appear as an excess at large 
E™ 1SS above the standard model expectation. Figure pi includes, as an example, the simulated 
distribution for a benchmark GGMwino model with a squark mass of 1700 GeV, a gluino mass 
of 720 GeV, and a total NLO cross section of 0.32 pb. 

For purposes of interpretation, we divide the data into six bins of Ej 1SS , indicated in Table [2 
For each bin. Table [2] lists the number of observed events, the number of predicted standarc 
model events, the acceptance for the benchmark signal model, and the number of background 
events introduced by the predicted signal contributions to the control regions, where this latter 
quantity is normalized to the corresponding signal yield. 

No significant excess of events is observed. An exclusion limit on the signal yield is derived 
at 95% confidence level (CL), using the CL S method Il55ti57l . For a given signal hypothesis, 
the six Ej' ss signal regions are combined in a multichannel counting experiment to derive an 
upper limit on the production cross section. The results, presented in Section [8} account for the 
possible contribution of signal events to the two control samples, which lowers the effective 
acceptance by 10-20% depending on the assumed SUSY mass values. 
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19.7 fb' 1 (8 TeV) 
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Figure 3: Distribution of £j ISS from the single-photon search in comparison to the standard 
model background prediction. The expectation from an example GGMwino signal model point 
is also shown. In the bottom panels, the ratio of the data to the prediction is shown. The 
representations of uncertainties are defined as in Fig. [2] 


Table 2: Observed numbers of events and standard model background predictions for the 
single-photon search. The signal yield and acceptance for the GGMwino model with ihj = 
1700 GeV and mg = 720 GeV, with a total signal cross section of c^lo = 0.32 pb, are also 
shown. The last line gives the additional number of background events, normalized to the 
signal yield, which is associated with signal contributions to the two control regions. 


E™ 1SS range (GeV) 

[100,120) 

[120,160) 

[160,200) 

[200,270) 

[270,350) 

[350,oo) 

Multijet 

991 ±164 

529 ±114 

180 ± 69 

96 ±45 

12 ±12 

9 ± 9 

ISR/FSR 

54 ±27 

73 ±36 

45 ±23 

40 ±20 

20 ±10 

15 ±7 

EW 

37 ±4 

43 ±5 

23 ±3 

19 ±2 

8 ± 1 

4 ± 1 

Background 

1082 ± 166 

644 ±119 

248 ± 73 

155 ± 50 

39 ±16 

28 ±12 

Data 

1286 

774 

232 

136 

46 

30 

Signal yield 

19 ±3 

53 ±5 

51 ±5 

82 ±7 

78 ±7 

67 ±6 

Signal acceptance [%] 
Background from 

0.3 

0.9 

0.8 

1.3 

1.2 

1.1 

signal relative to the 
signal yield [%] 

2.1 

5.0 

5.6 

9.9 

26.7 

13.5 


6 Double-photon search 

Events considered for the double-photon search are collected using triggers developed for the 
discovery of the Higgs boson in diphoton events Il58l - j60l . These triggers use complementary 
kinematic selections: 

• two photons with pj > 18 GeV, where the highest pj photon is required to have 
pj > 26 GeV, while the diphoton invariant mass is required to be larger than 70 GeV. 

• two photons with pj > 22 GeV, where the highest pj photon is required to have 
pj > 36 GeV. 
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In addition, each photon must satisfy at least one of two requirements: a high value of the 
shower shape variable Rg Il52l or loose calorimetric identification. For the targeted signals, the 
combination of the two triggers is found to be 99% efficient. 

In the subsequent analysis, at least two photon candidates with pj > 22GeV and |?/| < 2.5 
are required. Events are selected if the highest pj photon has pj > 30 GeV. Jets must have 
pj > 40 GeV and \rj\ < 2.5, with each jet required to lie a distance A R > 0.5 from an identified 
photon. Only events with at least one selected jet are considered. 

The background is dominated by multijet events, which mostly consist of events with at least 
one genuine photon. Due to the requirement of two photons in the event, the EW and ISR/FSR 
backgrounds are negligible. 

The razor variables Mr and R 2 lTl8l fl9l are used to distinguish a potential signal from back¬ 
ground. To evaluate these variables, the selected jets and photons are grouped into two exclu¬ 
sive groups, referred to as "megajets" HT9l . The four-momentum of a megajet is computed as 
the vector sum of the four-momenta of its constituents. Among all possible megajet pairs in an 
event, we select the pair with the smallest sum of squared invariant masses of the megajets. Al¬ 
though not explicitly required, the two photons are associated with different megajets in more 
than 80% of the selected signal events. 

The variable Mr is defined as 

Mr = \Z(Ip j1 I + Ip J 2 |) - (pi 1 + p'z) , (i) 


where p ] ' and p\ are, respectively, the momentum of the zth megajet and the magnitude of its 
component along the beam axis. The pj imbalance in the event is quantified by the variable 
My, defined as 


My = 


\ 


E™ iss ( 

IPt 

+ IptI 

) - p miss . 

[Pt + Pt) 


( 2 ) 


where pj is the transverse component of p ] ‘. The razor ratio R is defined as 


R = 


M r 

Mr' 


( 3 ) 


For squark pair production in R-parity conserving models in which both squarks decay to a 
quark and LSP, the Mr distribution peaks at Ma = (m~ — ?zz 2 SP )/?/Zq, where ??Zq (?«lsp) is the 
squark (LSP) mass. Figure [4] demonstrates that Mr also peaks for gluino pair production (left) 
and in the GGMbino model (right). 

The (Mr, R 2 ) plane is divided into two regions: (i) a signal region with Mr > 600 GeV and R 2 > 
0.02, and (ii) a control region with Mr > 600 GeV and 0.01 < R 2 < 0.02. The control region is 
defined such that any potential signal contribution to the control region is less than 10% of 
the expected number of signal events, producing a negligible bias on the background shape 
determination, corresponding to less than a 2% shift in the predicted number of background 
events for 20 expected signal events. 

The background shape is determined through a maximum likelihood fit of the Mr distribution 
in the data control region, using the empirical template function 

P(Mr) cx e~ fc ( MR “ M! 0", 


( 4 ) 
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Figure 4: Distribution of Mr in the double-photon search for the background model, derived 
from a fit in the data control region, and for the T5gg (left) and GGMbino (right) signal models. 
The background model is normalized to the number of events in the signal region. The signal 
models are normalized to the expected signal yields. 

with fitted parameters k, Mj*,, and n. The best-fit shape is used to describe the Mr background 
distribution in the signal region, fixing the overall normalization to the observed yield in the 
signal region. This implicitly assumes a negligible contribution of signal events to the overall 
normalization. We have studied the impact of the resulting bias and found it to be negligible 
for the expected signal distributions and magnitudes. The covariance matrix derived from the 
fit in the control region is used to sample an ensemble of alternative Mr background shapes. 
For each bin of the Mr distribution, a probability distribution for the yield is derived using 
pseudoexperiments. The uncertainty in each bin is defined by requiring 68% of the pseudoex¬ 
periments to be contained within the uncertainty band. 

This background prediction method is tested by applying it to a control sample of events in 
which jets are misidentified as photons, obtained by selecting photon candidates that fail the 
requirement on the cluster shape or the photon isolation. The remainder of the photon-selection 
criteria are the same as for the signal sample. In Fig.[5]we show the fit result in the control region 
(left) and the extrapolation to the signal region (right). 

The contribution of the EW and ISR/FSR backgrounds, characterized by genuine Ej' ss , is eval¬ 
uated from simulated events and is found to be negligible compared to the systematic uncer¬ 
tainty associated with the multijet background method, and is accordingly ignored. 

A signal originating from heavy squarks or gluinos would result in a wide peak in the Mr 
distribution. This is shown in Fig. [6j where a GGMbino signal sample is added to the control 
sample of jets misreconstructed as photons, and the background prediction method is applied. 
The contribution of signal events to the control region is negligible and does not alter the back¬ 
ground shape of Fig. [5] (left). The signal is visible as a peak at around 2TeV. 

Figure [7] (left) shows the result of the fit and the associated uncertainty band, compared to the 
data in the control region. The fit result is then used to derive the background prediction in the 
signal region. The comparison of the prediction to the observed data distribution is shown in 
Fig. [7] (right). No evidence for a signal is found. The largest positive and negative deviations 
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7 Signal model systematic uncertainties 
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Figure 5: Distribution of Mr in the double-photon search for a control sample of jets misre- 
constructed as photons (see text) in the control (left) and signal (right) regions. The data are 
compared to the 68% range obtained from a fit in the control region and extrapolated to the 
signal region (blue bands). The open dots represent the center of the 68% range. The rightmost 
bin in each plot contains zero data entries. The bottom panel of each figure gives the z-scores 
(number of Gaussian standard deviations) comparing the filled dots to the band. The filled 
band shows the position of the 68% window with respect to the expected value. 


from the predictions are observed for Mr > 2.3 TeV and 1.1 < Mr < 1.9 TeV, respectively, each 
corresponding to a local significance of ~ 1.5 standard deviations. 


7 Signal model systematic uncertainties 

Systematic uncertainties in the description of the signals are listed in Table [3] Differences be¬ 
tween the simulation and data for the photon reconstruction, identification, and isolation effi¬ 
ciencies are listed as Data/MC photon scale factors. The uncertainty associated with the parton 
distribution functions (PDF) is estimated using the difference in the acceptance when different 
sets of PDFs are used I6T1465I . Similarly, different sets of PDFs and different choices for the 
renormalization scales yield different predictions for the expected production cross section. 

Table 3: The systematic uncertainties associated with signal model yields. For the double¬ 
photon razor analysis, the contributions labeled as "shape" have different sizes, depending on 
Mr. 


Systematic uncertainty 

Single photon [%] Double photon [%] 

Data/MC photon scale factors 

1 

1-2 

Trigger efficiency 

2 

1 

Integrated luminosity [66] 

2.6 

2.6 

Jet energy scale corrections [67] 

2-3 

shape (bin by bin) 2-5 

Initial-state radiation 

3-5 

<1 

Acceptance due to PDF 

1-3 

1-3 

Signal yield due to PDF and scales 

5-20 

1-50 
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Figure 6: Distribution of Mr in the double-photon search for a control sample of jets misre- 
constructed as photons to which a simulated sample of GGMbino events has been added. The 
squark and gluino masses are respectively set to = 1400 GeV and nig = 1820 GeV, and the 
production cross section is fixed to a = 2.7 fb. The signal contribution is shown by the red his¬ 
togram. The representations of the uncertainty bands, data points, and the information shown 
in the bottom panel are the same as in Fig. [5] 

8 Interpretation of the results 

The result of the single-photon analysis is used to extract a limit on the production cross sec¬ 
tions of the GGM and SMS models. Comparing the excluded cross section to the corresponding 
predicted value, a mass limit is derived in the squark versus gluino mass plane. This procedure 
allows comparisons with previous results |[23l . In the SMS, the limits are derived in the gluino 
versus gaugino mass plane. The resulting cross section upper limits and the corresponding 
exclusion contours are shown in Fig. [8] 

Figure [9] shows the excluded mass regions and the cross section upper limits for the GGMbino 
and T5gg models obtained from the double-photon analysis. 

The single- and double-photon analyses are complementary with respect to the event selection 
and the search strategy. While the former is a multichannel counting experiment based on the 
absolute prediction of the standard model backgrounds, the latter uses kinematic information 
about the razor variables to perform a shape analysis. The best individual sensitivity is in the 
wino- and the bino-like neutralino mixing scenario, respectively. The double-photon analysis 
performs slightly better compared to the single-photon search in the bino scenario, because of 
the high-Hj trigger requirement in the single-photon selection. 


CMS 


_ 19.71b' 1 (8 TeV) 

Razor yy, > 1 jet 

Control sample 
+ Signal simulation (high R 2 ) 
Estimation in high R 2 ± o, 


GGMbino 




m~=1400 GeV 

q 

m~=1820 GeV 


shape 


9 Summary 

Two searches for gauge-mediated supersymmetry are presented: a search based on events with 
at least one photon and at least two jets, and a search based on events with at least two photons 
and at least one jet. The single-photon search characterizes a potential signal as an excess in the 
tail of the E™ ISS spectrum beyond 100GeV, while the double-photon search exploits the razor 
variables Mr and R 2 . These searches are based on pp collision data collected with the CMS ex- 
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Figure 7: Distribution of Mr for the control (left) and signal (right) regions. The representations 
of the uncertainty bands, data points, and the information shown in the bottom panel are the 
same as in Fig. [5] 


periment at a center-of-mass energy of \/s = 8 TeV, corresponding to an integrated luminosity 
of 19.7 fb \ No evidence for supersymmetry production is found, and 95% CL upper limits 
are set on the production cross sections, in the context of simplified models of gauge-mediated 
supersymmetry breaking and general gauge-mediation (GGM) models. Lower limits from the 
double-photon razor analysis range beyond 1.3 TeV for the gluino mass and beyond 1.5 TeV for 
the squark mass for bino-like neutralino mixings in the studied GGM phase space, extending 
previous limits H23l by up to 300 and 500 GeV, respectively. The limits from the single-photon 
analysis for wino-like neutralino mixings range beyond 0.8 TeV for the gluino mass and 1 TeV 
for the squark mass in the same GGM phase space, extending previous limits by about 100 and 
200 GeV. Within the discussed supersymmetry scenarios, these results represent the current 
most stringent limits. 
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Figure 8: Single-photon analysis 95% CL observed upper limits on the signal cross section and 
exclusion contours in the gluino-squark (top) and gaugino-gluino (bottom) mass plane for the 
GGMbino (top left), GGMwino (top right), T5gg (bottom left), and T5wg (bottom right) models. 
The thick red dashed (black solid) line shows the expected (observed) limit. The thin dashed 
line and band show the 68% CL range about the expected limit. The solid line quantifies the 
impact of the theoretical uncertainty in the cross section on the observed limit. The color scale 
shows the excluded cross section for each set of mass values. 
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